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^T) \ Abstract 

Let m > 1 be an arbitrary fixed integer and let Nm{x) count the number 
of odd integers u < x such that the order of 2 modulo u is not divisible by 
^ , m. In case m is prime estimates for Nm{x) were given by Miiller that were 

subsequently sharpened into an asymptotic estimate by the present author. 
' Miiller on his turn extended the author's result to the case where m is a 

. prime power and gave bounds in the case m is not a prime power. Here 

an asymptotic for Nfn{x) is derived that is valid for all integers m. This 
asymptotic would easily have followed from Miiller's approach were it not 
. for the fact that a certain Diophantine equation has non-trivial solutions. 

^ I All solutions of this equation are determined. We also generalize to other 

base numbers than 2. For a very sparse set of these numbers Miiller's 



l/^ • approach does work. 

o 

O '. 1 Introduction 

+^ ■ Let u be odd. Denote by l{u) the smallest natural number such that 

B '■ 2'(") = l(mod u). 

• r-j . The number l{u) is called the order of the congruence class 2(mod u). Let Njn{x) 

rS I denote the number of odd integers u < x such that m f l{u). It was shown by 

I Franco and Pomerance P[ Theorem 5] in the context of a generalization of the 

celebrated "3x + 1" -problem, that almost all integers u have the property that 

m\l{u), i.e. they established that, as x tends to infinity, 

Nm{x) = 0{x). 

The object of this note is to derive an asymptotic formula for Nm{x). Partial 
progress towards this goal was made by H. Miiller in his papers jElIZI- Let q > 2 
be a prime. Miiller [Oj showed that 

X -71 -r / \ 

< NJx) < 



X 



This was improved by the present author in [3], where he showed that 



^9/(9^-1) a; \^ V log a; 
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with Cg > a positive real constant. On his turn Miiller [Zj improved on this, by 
estabhshing the following generalization of the asymptotic estimate (^. 

Theorem 1 (H. Miiller). Let q > 2 be a prime and n > 1 fixed. Then 
N (r)-r ^ (l I r>^(^"g^"g"^'' 



log5^-"/{.^-i)^ V V log a; 
with Cqn a positive real constant. 

Miiller did not obtain an asymptotic result in case m is not a prime power. If 
m = p^ ■ ■ ■p^'', he notes that 

iVpji {x) < NUx) < N^i^ (x) + ■ ■ ■ + iVp=. (x) (2) 

and from this infers, in case m is odd, that 

OC -71 -r / \ 

< Nm{x) < 



log X \og^ X 

where a := max{C(pj, ej)|l < j < r}, resp. /3 := m.m{C{pj, ej)\l < j < r}, 
where for any prime q and integer n > 1 we define C{q, n) := q^~^ j {q^ — 1). (In 
fact, Miiller erroneously swaps 'min' and 'max' in his definitions of a and (3.) 

Note that the stronger result, in case m is odd, that xlog~^x -C N^^x) ^ 
xlog~'^x actually follows from Q. More can be said however: 

Lemma 1 Let m > 1 be odd. Suppose that there is only one integer j such that 
C{pj,ej) = P, then we have that Nm{x) ~ N^'^j^x) as x tends to infinity and the 

asymptotic behaviour of Nm{x) is given by Theorem^ 

Proof. W.l.o.g. we can assume that C(pi,ei) = (3. Now if there is no further 
j such that that C{pj,ej) = P, then by Theorem [T] each of the terms N^<h,{x) 
with > 2 is asymptotically of smaller growth than A^^^^i (x) and so the lemma is 
proved. □ 

Note that if there is more than one integer j such that C{pj,ej) = (3, then 
from together with Theorem^ the asymptotic behaviour of Nm{x) cannot be 
determined. We are thus led to the problem of determining pairs of integers j and 
k and primes pj and pk such that C{pj, Cj) = C{pk,ek)- The proof of the following 
result uses some arguments kindly provided by J.-H. Evertse and Y. Bilu. It 
follows that C(2,5) = C(5, 2) = C(3, 3) = 1/24, C(2,6) = C(7, 2) = 1/48 and 
C(5,3) = C(ll,2) = 1/120. 

Proposition 1 There are only finitely many solutions {p,a,q,(]) to 

with a and j3 integers and p and q primes with p < q. These solutions are: 
(2,5,3,3), (2,5,5,2), (2,6,7,2), (3,3,5,2) and (5,3,11,2). 
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Proof. It is easy to see that if there is to be a solution, we must have a > 2 and 
13 > 2. Assume w.l.o.g. that q > p. If (3 > 3, then q must divide either p — 1 or 
p + I. Since q > p it follows that q = p + 1 and so p = 2 and q = 3. This gives 
rise to precisely one solution: (2,5,3,3). So we may assume that [3 = 2 and we 
are reduced to finding the solutions of 

p\p' - 1) = - 1, (4) 

with 6 > 0. 

First assume that p = 2. Then we have to solve the equation — 1 = 3 ■ 2^*. 
As one of g it 1 cannot be divisible by 4, it follows that either g — 1g{1,2,3,6} 
or g + 1 G {1,2,3,6}. This gives rise to the solutions (2,5,5,2) and (2,6,7,2) 
(and no more). Thus we may assume that p > 2. There are two cases to be 
considered: 

1) p^\q + 1. Then we can write g + 1 = p^r, q — 1 = s (say) and so rs = p^ — 1. 
Thus p^r — s = 2 and hence r{p''r — 2) = — 1. Note that if 6 > 1, then we must 
have r > 1 and this gives rise to a contradiction. So 5 = 1 and hence we must 
have r{pr — 2) = p^ — 1. From this we infer that < r < p and, furthermore 
that — 2r = — l(mod p). It follows that r = [p + l)/2. Substituting this back we 
obtain — 2p — 3 = which only gives rise to the solution (3, 3, 5, 2). 

2) p^\q — 1. Now we have that g + 1 divides p^ — 1 and we obtain p^ < q — 1 < 
g+ 1 < — 1, which is impossible when 6 > 1. We can write g — 1 = pr, g+ 1 = s 
(say) and so rs = — 1. Thus s — pr = 2 and hence r{pr + 2) = p"^ — 1. From 
this we infer that < r < p and, furthermore that — 2r = l(mod p). It follows 
that r = {p—l)/2. Substituting this back we see that p^ — 6p + 5 = which only 
gives rise to the solution (5, 3, 11, 2). 

On collecting all solutions found along the way, the proof is completed. □ 

Since the Diophantine equation Q has non-trivial solutions we are blocked for 
certain integers m (for example for m = 2^ ■ 5^ = 800), in proving an asymptotic 
for Nm{x) by using (j21). We will provide an alternative to Lemma El and 
use this to obtain a more precise estimate for Nm{x) than can be provided by 
Miiller's method, and works for every integer m. At the same time we generalize 
to other base numbers than 2. 

2 Generalization to other base numbers 

First we will generalize the main result mentioned sofar. Theorem 1, to the case 
where the base number g is rational and not in {—1,0, 1} (an assumption on g 
maintained throughout this paper). We let uj{n) denote the number of distinct 
prime divisors of n and by Vp{n) denote the exponent of p in the prime factoriza- 
tion of n. Let S{g) be the set of integers composed only of primes p that do not 
occur in the prime factorization of g (i.e. of primes p such that Vpig) = 0). For 
each integer u in this set the order of g modulo u, ordg(u), is well-defined. We let 
Pg{d){x) denote the number of primes p < x with p G S{g) such that d\oTdg{p). It 
turns out that the set Pg{d) of primes p G S{g) such that d\ordg{p) has a natural 
density, which will be denoted by 6g{d). This density was first determined by 
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Wiertelak who derived a rather comphcated exphcit formula for it, a formula 
which was subsequently simplified by Pappalardi fH^ and streamlined further by 
Moree |3]. It turns out that 6g{d) is always a positive rational number. Part 1 
of the following result is due to Wiertelak part 2 to Moree As usual the 
logarithmic integral is denoted by Li(a;). 

Theorem 2 (Wiertelak P, Moree gl). 

1) We have Pg{d){x) = 5g{d)U{x) + O (^j-J-(loglogx)-W+3) . 

2) Under GRH we have Pg{d){x) = 6g{d)U{x) + 0(v^log'^('^)+^ x). 

In order to explicitly evaluate Sg{d), some further notation is needed. Write 
g = ±(70 5 where go is positive and not an exact power of a rational and h as large 
as possible. Let D{go) denote the discriminant of the field Q{^/go)■ (This nota- 
tion will reappear several times in the sequel.) Given an integer d, we denote by 
d°° the supernatural number (sometimes called Steinitz number), lipid Note 

that gcd(f,rf°°) = ripidP''"^''^- 

Definition. Let d be even and let eg{d) be defined as in Table 1 with 7 = 
max{0, U2{D{go)/dh)}. 

Table 1: eg{d) 



9\l 


7 = 


7 = 1 


7 = 2 


g>0 


-1/2 


1/4 


1/16 


g<0 


1/4 


-1/2 


1/16 



Note that 7 < 2. Also note that eg{d) = (-1/2)^^ if 50 
Theorem 3 (Explicit evaluation of 6g{d)). We have 



6g{d) 



d gcd(/i, d° 



n 

P\d 



e[id) 



(1 

1 + 3(1- 
1 + 3(1- 
1 

ll + %|(rf) 



sgn(^))(2'^2W 
sgn((7))(2'^^W 



with 



tf2\d; 

if2\ \d and D{go) \Ad; 



l)/4 

l)/A + eg{d) if 2\\d and D{go)\Ad; 

zfA\d, D{go)\M; 
%fA\d, D{go)\Ad. 



By N{x] g, m) we denote the number of integers u < x such that u G S{g) and 
m\oTdg{u). Note that N{x;2,m) = Nm{x). A straightforward generalization of 
Theorem 1 yields the following result. 



Theorem 4 Let n > 1 be fixed. Then 
N{x;g,q'') = Cgn{g) 



X 



O 



(log log x) 
logx 
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It is an easy consequence of Theorem |H1 that Sg{q^) = C{q,n) for almost all 
integers g. Furthermore, by Theorem IHl we find that 



17/24 if g = 2 and n = 1 

5/12 ifg = 2andn = 2 

2i-"/3 if g = 2 and n > 3 

C{q,n) otherwise. 



This evaluation allows us to formulate Theorem 1 for every prime power (and 
thus also for every power of two). It then follows from (2) that x\og~^^^"^^ x <^ 
Nm{x) <^ x\og~"'^^'^^ X, where 

7g(m) := mm{6g{p''/) \ 1 <j <r}. 

Lemma 1 can then be extended to all natural numbers m as follows: 

Lemma 2 Suppose that there is only one integer j such that 52{p^/) = 72(^), 
then we have that Nm{x) ~ N^^<^j{x) as x —>■ oo and the asymptotic behaviour of 

Nm{x) follows from Theorem^ 

We leave it as an exercise to the reader to show that the density of integers m 
satisfying such that 52 (p^^) = 72 1?^) for only one integer j exists and equals 

147497571941 



147916692000 



0.9971665 • • • 



An integer g such that there is only one integer j such that Sgijfj') = '~ig{rn) for 
every natural number m, we define to be a Miiller number. Miiller's inequality 
(2) generalizes of course to 

N{x-g,pl')<N{x-g,m)<N{x-g,pl') + --- + N{x]g,pl^). (5) 

The usefulness of Miiller numbers is apparent from the following result. 

Proposition 2 If g is a Miiller number and 5g{j)^-') = •yg^m), then we have 
Nm{x) ~ N^'^j{x) as X ^ 00 and the asymptotic behaviour of Nm{x) is given by 

Theorem |^ 

In the next section it will be seen that, unfortunately, Miiller numbers are very 
sparse. 



3 Some further Diophantine considerations 

The above discussion shows that if for (7 = 2 we want to cover odd integers m 
as well, rather than asking for the solutions of C{pi,ei) = C{pj,ej), we should 
be asking for solutions of S2{Pi') = 62{p^/), and indeed, more generally, for the 
non-trivial solutions (pi, ei,p2, 62) with pi 7^ p2 primes and ei, 62 > 1 of 

W)=W)- (6) 

A variation of the proof of Proposition 1 gives that in case g = 2 we only have 
the non-trivial equalities 52(2^) = ^2(3^) = 52(5^), 62(2^) = 62(7^) and ^2(5^) = 
52(11^). The more general situation is described by the following result. 
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Theorem 5 There are only finitely many quadruples {pi, ei,p2, 62) with pi < p2 
primes and 61,62 > 1 such that Sg{pl^) = Sg^p^^). These quadruples are given as 
follows (and there are no further ones): 

' (2, 5 - Ti - U2ih), 3, 3 - usih)) if 2^'^^ \ h and 3^ \ h 

(2,5-ri-z/2(/i),5,2-z/5(/i)) if2^-^^ \h and 5"^ \ h 

(2, 6 - ri - z/2(/i), 7, 2 - ujlh)) if 2^-^^ \ h and 7^ \ h 

(3,3-z/3(/i),5,2-z/5(/;.)) if 3^ ] h and5^]h; 

(5,3-z/3(/i),ll,2-z/ii(/i)) if5^ \h and 11^ \ h, 



where 



nig) 



^fD{go) = 
otherwise, 



and T2{g) = ri{g) + 



sgn(^) 



The associated values of 6g are 1/24, 1/24, 1/48, 1/24, respectively 1/120. 
Corollary 1 Let 

5i = {8000, 165373, 193600, 196000, 209088, 4002075, 4743200, 5122656}. 

Let S2, respectively 5*4 be Si, but with all even numbers in this set divided by 2, 
respectively 4. For a given number g equation (jHl) has only trivial solutions, that 
is g is a Miiller number, if and only if h is divisible by a number from S2T2(a) ■ 



Example. The numbers 
bers. The number g = 
Theorem El 



22000^ 2^0°° and 3^°°° are examples of small Miiller num- 

^ W by 



^4000 jg j-^Q^ ^ Miiller number, since Sg{2^ 



For the convenience of the reader we include a table where ti 
T2 = T2{g) are given. 

Table 2: Ti{g) and T2{g) 



Ti{g) and 





D{go) 7^ 8 


D{g^) = 8 




ri 


T2 




^2 


g>Q 








1 


1 


g<Q 





1 


1 


2 



Clearly Miiller numbers are very sparse. Indeed it is not difficult to quantify 
this: there are positive constants 64. and c_ such that the number of positive 
Miiller numbers up to x grows as c+x^/^^"" and the number of negative Miiller 
numbers with absolute value not exceeding x as cx^^'^^^^ as x 00. 

The proof of Theorem El is unfortunately not very instructive as many cases 
have to be distinguished. As its level of difficulty is comparable to that of the 
proof of Proposition 1, we leave it to the interested reader. 

In order to derive Corollary 1 we notice that the set of those integers h such 
that none of the 5 solutions in Theorem occur, are the natural numbers of the 
form {hi) U ■ ■ ■ U (/is), where {hi) denotes the ideal generated by hi and hi, . . . ,hs 
are certain integers. A priori the 5 solutions yield 2^ = 32 potential generators 
(for example lcm(3^, 5^, 7^, 5^, 11^) = 4002075), of which some turn out to be 
divisors of others (the latter ones can thus be left out from S'graCg)). 
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4 A more refined estimate for N(x] m) 



We sharpen Miiller's estimate (0) to an equahty (given in Lemma Ej) and use this 
to obtain a more refined estimate for N{x;g,m) (and thus Nm{x)). As before 
we let m = pi'^^ ■ ■ ■pr^'' denote the factorization of m. By K,{m) = pi ■ ■ ■ p^ we 
denote the squarefree kernel of m. A divisor (i of m is said to be unitary if 
gcd{d,m/d) = 1. Note that the non-trivial unitary divisors of m come out as 
'blocks' from its factorization. Let N'{x;g,m) denote the number of integers 
u < X from S{g) such that p^^ \ ordg(M) for 1 < j < r. We let P'{x;g,m) be 
similarly defined, but with the phrase 'integers m < x' replaced by 'primes p < x\ 
Note that N'{x; g,pf) = N{x; g,p''/). 

Lemma 3 We have 

N{x; 5(, m) = - ^ /i(K(d))iV'(x; g, d), 

d>l 

where by d^m we indicate that the sum is over the unitary divisors d of m. 

Proof. If M ^ S{g)^ the integer u is not counted on either side of the purported 
identity. If m < x is in S{g) we will show that it is counted with multiplicity one 
on both sides of the purported identity. If m|ordg(ti), then u is counted neither 
on the left nor on the right hand side. If m f ordg(M), we can assume w.l.o.g. 
that pj^^ \ oidgiu) for 1 < j < /c and that pl^loidgiu) for s > k. By definition u 
is counted with multiplicity one on the left hand side. The contribution of u to 
the right hand side is 

- ^ K<d)) = - J2 Mrf) = Mi) = i, 

d>l d>l 

where we used the well-known identity J2d\n^^id) ~ 0' which holds for every in- 
teger n > 1. □ 

Example. We have N{x; 2, 12) = N^ix) = N^ix) + N^{x) - N'{x; 2, 12). 

Put5^(d) = E,#„^(K(j))53(j)- 
Lemma 4 We have 

1) We have P'{x-g,d) = 5;(rf)Li(x) + O (^(log logx)-W+3y 

2) Under GRH we have P'{x; g, d) = 5'g{d)U{x) + 0(v^log'^('^^+^ x). 

Proof. By inclusion and exclusion (or an argument as in the proof of Lemma E)) 
we see that P'{x;g,d) = /i(K(j))Pg(j)(x). Then use TheoremEl □ 

j#d 

Remark. Wiertelak studied P'{x] g, d) in depth in case d is squarefree. Note 
(as did Wiertelak) that if d is squarefree, then P'{x] g, d) counts the number of 
primes p < x with I'p^g) = such that the congruence g'^''^ = g{mod p) has a so- 
lution y. In general P'(x; g, d) counts the number of primes p < x with i^pig) = 
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such that the congruence g'^^ = g^{modp) has a solution {r,y) with r dividing 
d/K,{d). 

The following result gives some insight in the properties of Sg{d). 

Proposition 3 Suppose that d > 1 and di > 1 are natural numbers. 

1) We have 6g{ddi) < 6g{d). If g > 0, then we have 6g{ddi) < ^6g{d). 

2) We have 

rSniddi) , , - , ,6a(ddi) , , , 5 

sup{ '' \di > 1} = 1 and sup{ ,\' \g > 0, di > 1} = -. 
6g{d) 6g{d) 6 

3) We have Sgi^pl^p^^) < Sg{pl^), with pi and p2 primes. 

4) We have Sg{pl^ ■ ■ ■p"/) = (1 - ■ ■ ■ (1 - ^givl'')), with pi,. . . ,pr distinct 
odd prime numbers. 

5) We have min{l - 5'g{j)\jifd, d>l} = -fg{d). 

Proof. 1) li g > 0, then one checks that e[{ddi) < je[{d) and from this one easily 
infers that 6g{ddi) < ^6g{d). The general case requires some case distinctions 
and is left to the interested reader. 

2) We have, e.g., limg^oo ^-52= (6)/5_52<! (3) = 1 and 53(6) = 553(3)/6. Now invoke 
part 1. 

3) On invoking part 1 we find that 6g{pl'pl^) - ) = 5g{pl'p2^) - 6g{p2^) < 0. 

4) Considered as a function of d, Sg{d) is a multiplicative function on the odd 
integers by Theorem IHl From this and the expression of 6g in terms of 6g, the 
result then follows. 

5) On noting that 6g{ddi) < Sg{d) and using part 3 we see that 

min{l - 6'g{j)\ji^d, J > 1} = min{l - 6'gip'r^)\l < j < r} = jg{d), 
where we used that 1 — Sg^p^') = dgijP^). □ 

Next we will compute N'{x;g,d) from P'{x;g,d). To this end we will make 
use of the fact that the integers counted by N'{x; g,d) as x ^ 00 can be related 
to completely multiplicative sets. Recall that a set S of natural numbers is 
said to be completely multiplicative if its characteristic function x is completely 
multiplicative, that is satisfies xi^u) = x{^)x{y) for all natural numbers x and 
y. Thus S is completely multiplicative when x-yESiSx&S and y E S. The 
following result is a simple generalization of Lemma 1 of Miiller [Z| and easily 
follows from the observation that ordg('Ui'U2)|gcd(-Ui, 'U2)lcm(ordg('Ui), ordg('U2)), 
cf. p. 

Lemma 5 Let m = pi^"^ ■■■pr'^''. The set of integers u from S{g) such that 
gcd{u,m) = 1 and p^-' \ oidg{u) for 1 < j < r is a completely multiplicative set. 

Let N"{x; g,m) denote the number of integers u < x from S{g) such that 
gcd(M,m) = 1 and p^^ \ oidgiu) for 1 < j < r. The integers thus counted (as 
X — 00) form a completely multiplicative set by Lemma El and can be counted 
using the following proposition. 



8 



Proposition 4 1) Let S be a completely multiplicative set of natural numbers 
such that for x ^ oo we have 



(7) 



holds, where r > and 7 > 0, with t real and 7 an integer. Then we have for 
X — > 00 that 



n<x 
neS 



where cg denotes a positive constant which only depends on the set S. 
2) // dZj) holds with error term 0(xlog~^~'''^ x) for some 71 > 0, then 

S{x)=x 6^1og"-^-''x + 0(slog"-^-^^+'x), 

0<i/<7i 

where &o(= c^), &i, . . . are constants. 

Proof. The first part is proved in Moree [2]. The second part is a special case of 
Theorem 6 of 0. □ 

Let m > 1 be fixed. Then, by the latter proposition and Lemmas 0] and 

N [x;g,m) = c^ig) — . . . 1 + 

^ ' ' log^-^o^"^^ X V V log a; 

where c'^{g) is a positive constant. 

As in |7] we conclude that there exists a finite set of numbers (actually prime 
powers) ni, . . . , (depending on g and m) such that 

s 

N'{x-g,m) = J2N"{-;g,m). (9) 
It thus follows on invoking the estimates (jH)) and © that 



X 



(log logx)'^(™)+^ 



A"(:.; 9, rn) = U9)^^^r^^ ^1+0 [^^^-^^ j j . (10) 

On noting that N'{x; g,p^j') = N{x] g,pj) and that 5'g{pJ) = 1 — ^g{pj)i we see 
that the estimate (fTIHl generalizes Theorem |3] 

Using (jlOj) we can now generalize Proposition 2. 

Proposition 5 Assume that m = pl^ ■ ■ -p^^ with 5g(Pi^) < ■ • ■ < 5g{p';^). Then 
we have N(x; g, m) ~ '^^"^^'^^ N(x; g,Pi'). It follows that in particular we have 
N{x; g, m) ~ ex log"*^^*-^! ^ x for some c > as x —>■ oo. 
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Proof. By Lemma El (fTO|) and part 5 of Proposition 01 we have N{x]g,m) ~ 
jyj=i ^ i^i 9 y pT) ■ Since by Theorem we cannot have that Sg{pl^) = Sgijh^) = 
^giPz'') — ^givT)) the result follows on invoking (fTUI) again. □ 

The main result of this paper makes the latter result more precise: 

Theorem 6 Let | j#m, j > 1} = . . . , 54 with 5i < 82 < ■ ■ ■ < 5s. 

1) Then 5i = 'Jgim) and 

where the leading (asymptotically dominating) term is the one with j = I and 
Ci, . . . ,Cs are constants with Ci > 0. 

2) Under GRH it is true that for each integer t > we have 

s t 

N(x;q,m)=} > Cj k m \- 0( rri — )• 

log'^+^^ x ^\og'+^x' 

j=l fc=0 ^ " 

The constants Cj and Cj^k will depend on g and m in general and so do the implied 
constants in the error terms. 

Proof. 1) Follows on combining Lemma El pO|) and part 5 of Proposition 01 

2) The proof is similar to that of part 1, except that instead of part 1 of Lemma 

0]we use part 2 and that instead of part 1 of Proposition 0] we use part 2. □ 
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